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Abstract 

Chiral symmetry is restored at high density, quarks become nearly massless and pion, 
the Goldstone of the symmetry breaking decouples from the quarks. What happens at 
high density is important for finding the density dependence of Strange Quark Matter 
(SQM), - which in turn is relevant for understanding the structure of compact stars. 
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1 Introduction 



We investigate what happens to /„-, at high densities. In our convention, /„. is defined as 
follows (vacuum value of f w ~ 93 MeV) : 

<0\A*(x)\7r\q)> = iq»5 ab U{q 2 )e-^. (1) 

The Hellmann-Feynman theorem, applied to a nuclear many body model, gives the quark 
condensate in nuclear matter at high density [1] . Coupling this with theoretical one of Nambu 
and Jona-Lasinio (NJL) [2], one can extract f^ins), where n B is the baryon number density. 
Chiral Symmetry breaking and pion properties was discussed in the framework of NJL model 
by Bernard [3]. There was follow up of the work on /^(tib), using the NJL model, by Bernard, 
Meissner and Zahed [4] and more recently by Caldas [5]. 

Low temperature QCD sumrule results also give upto % ~ 4 no where no is normal 

density [6]. 

Again density dependent quark masses, used for SQM calculations [7], can be used to fix 
the parameters of the NJL model. This in turn enables one to get the pion coupling to the 
QCD vacuum f n (n B ). 

The quark mass is given in the SQM [7] as : 

M* = m i + M Q sech (^-) , i = u,d,s. (2) 

where n B = (n u +nd+n s )/3 is the baryon number density; no = 0.17 fm~ 3 is the normal nuclear 
matter density; n u , n<i, n s are number densities of u, d and s quarks respectively and N is a 
parameter. The current quark masses (m,) are taken as : m u = 4 MeV, m d = 7 MeV, m s = 
150 MeV. Mq is the constituent quark mass taken around ~ 325 MeV according to latest 
version of the model [8]. 



2 Nuclear matter model 

In the relativistic a — u models of nucleon matter it is found that the quark condensate can be 
estimated using the Hellmann-Feynman theorem and this was investigated in detail in [9, 1]. 
Interestingly, the title of [1] also referred to a decoupling, - that of the nucleon mass and the 
quark condensate in the medium. The Walecka model, the pioneering one, implies an effective 
quark condensate that increases with density. This is contrary to common belief. The newer 
Zimanyi-Moskowski (ZM) model, has an edge over the Walecka model in satisfying the criterion 
that the quark condensate falls with increase in density as shown in [1]. 

Further relevance of the ZM has been recently pointed out by Sinha et al., who have shown 
that the velocity and the incompressibility of the ZM model also match onto a quark model 
[10]. 

According to Hellmann-Feynman theorem [11, 12, 9] 
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< V(A) | ±H(\) | V(A) > = A < ^(A) | //(A) | V(A) >, (3) 

where H(\) is any hermitian operator depends on a real parameter A and |^(A) > is a normal- 
ized eigen vector of H(X). 

In QCD the Hamiltonian density is given by 

7-Lqcd = H + 2m q qq, (4) 

with the major part being the chirally symmetric Ho. Here m q is quark mass and q is the quark 
field. 

Making the identification H — > / d 3 x Hqcd and A — > m q one finds the Hellmann-Feynman 
theorem as: 

2m g < ^(A) | y rf 3 x gg | ^(A) > = ™<?^- < ^(A) | ^ d 3 a: Hqcd I ^(A) > . (5) 

The above equation may be applied to nuclear matter and vacuum with | ip(\) > = \ n B > 
and | ip(X) > = | vac > respectively. Here | % > denotes ground state of nuclear matter 
at rest with nucleon density tib and | vac > denotes the vacuum state, taking the difference of 
the above two cases and keeping in mind the uniformity of the system, one gets 

2m q (< qq > nB - < qq > vac ) = m q j^- (6) 

where ub is the number density in nuclear matter. 

Here in general < Q > nB = < n B \ | % > and < Q > vac = < vac \ Q \ vac > 
notations have been used for an arbitrary operator Q. 

The energy density £ of nuclear matter is given by 

S = n B M N + 58, (7) 

where 88 is the contribution to energy density from the nucleon kinetic energy and nucleon- 
nucleon interaction energy. 88 is of higher order in the nucleon density and is empirically small 
at low density. 

At low density the quark condensate can be related to the nucleon a term a N which may 
be defined as [13] 

°n = \ E(< N I [Q% [Q%H QCD ]] \N>-<vac\ [Q% [Q a A ,H QCD ]] | vac >), (8) 

6 a=l 

where Q a A is axial charge, Hqcd QCD Hamiltonian and | iV > is state vector of nucleon at rest. 
Alternatively ctn can be expressed as: 
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(Tat = 2m q J d 3 x (< N \ qq \ N > — < vac \ qq \ vac >). (9) 

where 

dM N 

a * = m 9-r— • ( 10 ) 



dm q 



Hence Eq. (6) can be written as (using Eq. (7)) 



. . dM N d 58 d 58 , 1lN 

2m q {<qq> nB - <qq> vac ) = m q n B — V m q - = n B a N + m q - . (11) 

dm q dm q dm q 

Assuming translational invariance which makes quark condensate constant one can define 

a A = 2m q V (< qq > riB - < qq > mc ). (12) 

Using Eq. (11) 



°A , . d 5{8/n B )\ , v 

^ = T = aN [ 1 + ^M^) (13) 

Now from Gell-Mann— Oakes— Renner relation we know 

2m q < qq > vac = - m 2 Jl (14) 
and f n being the pion mass and pion decay constant respectively. From Eq. (6) 

^^ = l-n B ^J- (15) 
< qq > 

vac 

In the ZM model the Lagrangian describes the motion of a baryon with an effective mass 
instead of bare mass. This information goes to modify the scalar coupling constant making it 
density dependent while the vector coupling remains the same. In contrast with the Walecka 
model < qq > nB /< qq > vac goes down with density [1]. 



3 The QCD sumrule method 

This is a very elegant method devised by Shifman, Vainshtein and Zakharov [14] and consists 
of equating the coupling of an interpolating Lorentz invariant current for a meson or a baryon 
- with proper spin, parity and isospin degrees of freedom - to quark- ant iquark for meson and 
three quark for baryons. The quarks or antiquarks are then allowed to mix into the QCD 
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vacuum - which have condensates of quark- ant iquark and gluons - and also exchange gluon 
lines through operator product expansion (OPE). Starting at high momentum transfer for 
finding the coefficients of the OPE by Borel transform one finds a 'window' where the sum 
rule becomes independent of he Borel mass parameter. The condensate values picked up from 
one set, say the p meson can be used for all the meson or baryon sumrules. For meson baryon 
coupling one has to go over to three point functions which is more complicated but straight 
forward in principle. Reviews are available by Reinders, Rubinstein and Yazaki [15] and Dey 
and Dey [16]. 

Working out the density and temperature dependence of the pion-nucleon coupling constant 
{9-kNn), within the framework of QCDSR techniques, Dey and Dey [6] deduced the f n to be 
about half its value (44 MeV compared to 81 MeV) at four times normal density 1 . This agrees 
with the estimate of the present paper using the NJL model. The sumrule model predicted 
that the Goldberger-Treiman relation g^NN = M N \^2/ f n is independent of density [17] and 
this was confirmed in a later calculation by [18]. 

4 Quark mass used in stellar calculation 

Early suggestion of a cosmic separation of phases of hadronic and strange matter led to inves- 
tigations properties of strange quark star, but were not very successful. This was because the 
star with maximum mass had a radius of about ~ 9 — 10 kms and this is comparable to that of 
a neutron star. One could not distinguish between the two. The density dependence of quark 
masses was not considered in these early models. At high density there is chiral symmetry 
restoration (CSR) and the masses approach the current quark mass values. 

By putting CSR, in a simple tree level large N c model [7] , one can set up an equation of state 
(EOS) and seek to explain the properties of compact stars Her X— 1 and 4U 1820—30. Li and 
others used this EOS to explain the properties of SAX J1808.4-3658 or 4U 1728-34 [19, 20]. 
Compact stars are assumed to be composed of (u, d, s) matter that is very dense (typically 4.6 
(surface) to 15 (core) times the normal nuclear density). In the model (u,d) matter has less 
binding per baryon E/A, compared to Fe 56 and (u, d, s) matter has more. 

It is interesting to note that many X-ray emitters are rotating and shows periodicity. Only 
recently, however, six sources were discovered starting with SAX J1808.4— 3658 (1998), which 
are accreting millisecond X-ray pulsars and an important question is raised by Wijnands [21] 
: why are those compact stars different from others for which no pulsations have been found? 
Perhaps, he comments, new ideas need to be explored to explain these six sources. Stability 
of the star may be a crucial point in resolving this issue, according to the present authors and 
the use of (u, d, s) matter with restored chiral symmetry may help. We must mention that 
the model leads to stars which are very stable as shown by Sharma et al. [22], by matching 
the external Schwarzchild metric to a realistic one at the boundary of the star. For details we 

1 f 7T was normalized to vacuum value 130 MeV in [6] and is readjusted here by the factor y/2. It is somewhat 
low in a nuclcon which already has a substantial hadron density. 
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refer the reader to [22]. Strange star models, with the above EOS, are also very stable when 
rotating fast, as shown by Gondek-Rosihska et al. [23] and Bombaci, Thampan and Datta [24]. 
The density dependence of the strong coupling constant a s in this model was explored using 
the simple Schwinger-Dyson expansion advocated by Bailin, Cleymans and Scadron in Ray et 
al. [25]. 

Further, there are other interesting applications of this model enumerated below : 

1. X-ray superbursts lasting for several hours thrice in 4U 1636—53 and once (so far) in KS 
1731—260 [26], and also the phenomenon in general, seen in 7 stars altogether. 

2. Occurrence of two quasi periodic peaks in the X-ray power spectrum model of 4U 1636—53 
and KS 1731-260 [20, 27] and other stars. 

3. Absorption in 1E1207.4-5209 [28] and emission [29] in various stars like 4U 0614+091, 
2S 0918-549, 4U 1543-624, 4U 1850-087 from surface compressional modes. 

In addition the interesting model of quark nova of Ouyed et al. [30], employs the idea 
of contraction of normal matter when it is converted to (u, d, s) matter of the above model. 
Gravitational energy from matter falling onto a compact core, formed during a supernova 
explosion and consequent generation of a core remnant, can lead to gamma ray after glow 
according to [30]. 

The density dependent quark mass used in the (u, d, s) matter is used to generate the pion 
coupling to quarks in the present paper. 



5 The Nambu Jona-Lasinio model 

We recall that in the model of Nambu and Jona-Lasinio, one can calculate the quark mass 
M*, /„., the quark condensate < qq > for a given coupling G, - following the equations below 
in terms of a cut off A of 631 MeV (see [31]) : 



M* = m + 4G{N c N f + \)W f-j^L^ (16) 



(2tt) 3 E 

f 2 = N r M* 2 f A d " P — (17) 
T * c J (2vr) 3 £ 3 1 ' 

- d^v I 

<qq> = <uu> = <dd> = — 6M* / 7—^- (18) 



(2tt) 3 £ 

Knowing the NJL coupling G, one can therefore relate the quantities M*, f n and qq. We 
assume that G varies with density and find it (1) by fitting it to /V in the QCDSR for which 
we do not need the NJL model, (2) using < qq > in the a — uj nuclear matter model and 
(3) by fitting density dependent (u,d,s) quark mass in equation (16). /,From (17) and density 
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Table 1: Variation of /„., G, and < qq > with density ratio (n = 0.17/ fm 3 ). 



n B /n 


u 

(MeV) 


G 

(MeV- 2 ) 


< qq > 1/3 
(MeV) 


1 


90.9227 


4.936 x 10~ 6 


-243 


2 


85.8209 


4.682 x 10~ 6 


-232 


3 


77.5264 


4.389 x 10~ 6 


-223 


4 


67.1299 


4.13 x 10~ 6 


-208 


5 


56.1471 


3.922 x 10~ 6 


-191 


6 


45.8124 


3.755 x 10- e 


-174 


7 


36.8022 


3.606 x 10~ 6 


-159 


8 


29.3312 


3.458 x 10~ 6 


-144 


9 


23.3414 


3.294 x 10~ 6 


-131 


10 


18.6502 


3.103 x 10~ 6 


-120 


11 


15.0386 


2.877 x 10~ 6 


-110 


12 


12.2944 


2.616 x 10~ 6 


-102 


13 


10.2314 


2.323 x 10 _e 


-95 


14 


8.69422 


2.01 x 10~ 6 


-89 


15 


7.55736 


1.692 x 10~ 6 


-84 


16 


6.72196 


1.386 x 10~ 6 


-80 


17 


6.11141 


1.107 x 10~ 6 


-78 


18 


5.66718 


8.645 x 10~ 7 


-75 
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dependence of G, f w and the corresponding quark condensates are obtained and are tabulated 
in table 1. 

At high density, nucleon mass decreases very much with /„. in the Skyrme and other models 
and the nuclear radius becomes so large that there is no point in talking of a 'confined' nucleons, 
the quarks are percolating. 

Fig. (1) shows that the a — uo model predicts a zero /„. at about ~ 4p - The QCDSR fall-off 
is also sharp compared to SQM. We can thus claim that the CR in SQM is mild. The full % 
dependence is shown in Fig. (2) where the density dependence of NJL coupling G is also shown. 
Our result checks with [4]. For example, for number density five times no the value of f n is 
about 60 MeV. A much more mild density dependence of f n is implied by Caldas [5] who display 
a number like 80 MeV. It will be very interesting to see if the photon width increase, predicted 
in this paper, is indeed found in heavy ion collisions. The photon momentum resolution of 
STAR experiment does not allow any decisive conclusion about the possible enhancement of 
the 7T° width, - for details see [5] . 

n i i i 

□ 

□ 

o □ 

+ 

o □ 

o □ 

+ 

o 



1 2 3 4 5 

n B /n 

Figure 1: Density f n from different models upto (~ 4p) : diamonds corresponds to QCDSR results, + 
corresponds to the nuclear matter model of ZM, squares correspond to the SQM. 

For future use we have fitted all the quantities by the equation 

y = a±x + (I2X 2 + a^x z + a^x A + a 5 x 5 + a 6 x 6 + a 7 x 7 + a$x 8 (19) 

where y represents the variables ( f n , G, < qq > respectively ) and x is the density ratio 
ub/tio. The coefficients for each quantity are tabulated in table 2. 

6 The problem of relativist ic heavy ion collisions (RHIC) 

Recently exciting new results have been reported by several groups from the gold on gold nu- 
clear collisions in Brookhaven. It appears that there is thermalization and a high temperature 
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2 4 6 8 10 12 14 16 



n B /n 

Figure 2: Predictions from density dependent quark mass of the SQM upto high density : density (p) depen- 
dence of f v (full line), p dependence of the NJL coupling constant G (dotted line ). 



Table 2: Coefficients for density expansion of /„-, G and < qq >. 



coef 


U 


G 


< qq > 


Cli 


164.51 


9.025 x 10" 6 


-2.65893 x 10 7 


a 2 


-104.31 


-5.832 x 10" 6 


1.72333 x 10 7 


a 3 


31.08 


1.780 x 10~ 6 


-5.11645 x 10 6 


0,4 


-5.21 


-3.006 x 10" 7 


849758.0 


a 5 


0.513 


2.959 x 10" 8 


-83014.7 


a 6 


-0.029 


-1.690 x 10" y 


4724.5 


a-j 


-0.0009 


5.187 x 10" 11 


-144.734 


a 8 


-1.156 x 10" 5 


-6.611 x 10~ 13 


1.843 
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is reached. The problem with the experimental results is that although the system is not de- 
scribable by hadronic models, the nearly non-interacting quark gluon model also does not seem 
to work. In the language of the protagonists ' the interpretation of current data relies heavily 
on theoretical input and modeling, in particular on the apparent necessity to include partonic 
degrees of freedom in order to arrive at a consistent description of many of the phenomenon 
observed in experimental data. Seen from a purely experimental point of view this situation 
is somewhat unsatisfying, but probably not unexpected, not avoidable, considering the com- 
plexity of the reaction and associated processes [32].' Quoting another group to conclude 'the 
data from RHIC collisions provide strong evidence for the creation of high energy density, low 
baryonic chemical potential, medium which can not simply be described in terms of hadrons 
and whose constituents experience significant interactions with each other [33].' 

In conclusion, from high temperature RHIC data, it is not clear that either of the features 
of QCD like chiral symmetry restoration (CSR) or asymptotic freedom (AF) is actually realized 
due to the complexity of the system and the system may display strong interacting coherent 
partonic interactions. The system that one can observe in stars may in fact yield a clearer 
signature of CSR and AF. We are grateful to the referee for allowing us to comment on this 
feature. 

In the next section we shall discuss the nature of the density dependence that one expects 
from heuristic considerations given by various authors. 



7 Discussion 

In the model [34], the radius of the pion is : 

R n = 0.4^ IU (20) 

where z is the probability of finding a purely qq component in the pion. The decrease of f n 
with increasing density signifies increase in the radii of the hadrons. This in fact ultimately 
leads to the percolation of the quarks. Assuming the nucleon radius Rjy = c MeV fm / f n , 
in [6] the constant c is adjusted to get the radius of the nucleon at normal density : 

R N = 86.12 MeV fm/f n . (21) 

One can review QCD scales following Bailin, Cleymans and Scadron [35]. 

m dyn ~ A Ms e^ ~ 300MeV (22) 

where the minimal subtraction overall energy scale of QCD, M^ s ~ 250 MeV for the 3 flavour 
case. This is close to 325 MeV of (2). One can go on to get 

U = ^rn dyn » 87MeV (23) 
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and the string tension 

a ~ ^m dyn « WOMeV. (24) 

As density increases /„- and quark condensate decreases with md yn and this is borne out by 
the NJL model of the present paper. 

f n is a parameter in chiral models, the pioneering one being the Skyrme model, 

C = \fl TridpUdJP) + (l/32e 2 ) Tr(d^UU\ d v UU^f (25) 

o 

where it is the only parameter depending on temperature and density [36]. The consequences 
of f w (p) was first analyzed by Rho [37] and Meissner [38]. 

The nucleon mass can be calculated using the Skyrme model : 

M N = Un/ey/2 (26) 

where \x = 73.0 is an integral over the chiral angle of the Skyrmion [39] and e is the dimen- 
sionless Skyrme parameter taken to be 5.78. 

In this context it is interesting to emphasize the suggestion by Dosch and Narison [40] , from 
QCD Sum Rule (QCDSR) method, that e is independent of the quark condensate. Based on 
this [36] found that indeed the parameter e, being independent of temperature and density, 
could in fact be 2tt, as suggested by Skyrme to represent a spin current. Incidentally the 
nucleon radius Rn is proportional to its inverse Rn = c MeV fm / f n where c is a constant. 
In general, all chiral model properties scale with f n , as in the Skyrme model. 

These values have the support of tentative observations made for compact stars. The 
importance of the results can be anticipated, since a convincing proof for the existence of such 
compact stars may soon emerge, from the copious flow of recent astrophysical observations. 

In particular it will be interesting to see if there is any change in Ouyed's model for Skyrmion 
Star [41] with a density dependent f w . 



8 Conclusions and summary 

We have calculated the variation of the pion coupling f n (p) with density in the Nambu 
Jona-Lasinio model and it is satisfying to see that this matches with expectations of other 
models. f n (p) and the constant G are parametrized as polynomials of density in the hope that 
the results may be used in future calculations. 

In summary, we have calculated the pion coupling constant /„. from the density dependent 
(u, d, s) masses employed in compact star models and the results are qualitatively matching 
with other models namely (1) QCD sumrule and (2) nuclear matter models. Results may be 
useful for chiral models where use of f n (p) will produce significant difference at high density. 

To conclude, in our opinion, observations on high density matter - perhaps possible in 
compact stars in an indirect manner - may yield signatures of asymptotically free and nearly 
chirally symmetric matter. These signatures are elusive in present day RHIC data. 
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